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Abstract 

We present a new integral transform called the Generalized Borel Transform (GBT) and show how 
to use it to compute some distribution functions used to describe the statistico-mechanical behavior 
of macromolecules. For this purpose, we choose the Random Flight Model (RFM) of macromolecules 
and show that the application of the GBT to this model leads to the exact expression of the polymer 
propagator (two-point correlation function) from which all the statistical properties of the model 
can be obtained. We also discuss the mathematical simplicity of the GBT and its applicability to 
polymers with other topologies. 
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I. INTRODUCTION 



Consider a group of n small molecules connected in a sequential manner such that each 
small molecule has only two nearest neighbours with whom it forms a chemical bond ex- 
cept for the small molecules at the end of the chain that have only one connected nearest 
neighbour. This chain is called a polymer. In this particular case, the chain ends do not 
form a chemical bond thus, the polymer is said to have a linear topology (linear polymer). 
If the ends were to form a chemical bond, then the polymer is said to be a ring (cyclic) 
polymer. Another well known topology is the star topology where m linear polymers grow 
from the same point and they are connected at this point. This defines an m-arm star 
polymer. Combinations of these three topologies or new topologies define more complex 
macromolecules. 

Another characteristic of polymers is that they can be made of two or more different 
kinds of small molecules. 

We conclude this section by speculating about other possible models for which the GBT 
might provide some useful insight. But, before we proceed, we should note that extensions of 
the present calculation to other polymer topologies like rings and stars are straightforward. 
For example, in the case of rings the cyclic topology of the polymer generates an extra 
constraint in the polymer propagator of the form 

which can be easily handled by the GBT. 

A generalization of Eq. (|23| that takes into account more complex systems like copoly- 
mers and the effects of external fields has the mathematical expression given by 



P(R,n,{p«})= /rf{R^n^Efe Q r;(R,))j5^R,-R^ exp (—u ({R,-})) , (2) 

where s is the total number of chemical species forming the copolymer, is the probability 
of finding the j-th segment in the a-th chemical species and uj ({R 3 -}) has the mathematical 
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form 



w ({R,}) = J> (R,-) , (3) 

77 (R) being any function. In particular, the effects of external vectorial (rj (R) = — F • R) 
and quadrupolar (r) (R) = QijRiRj) fields can be studied exactly. In general, the statistical 
properties of the model defined by Eq. (J2J) can be computed exactly using GBT. 

In recent years, new experimental tools have been developed to study the behavior of sin- 
gle macromolecules. These tools have given researchers the ability to manipulate polymers 
individually. Perhaps, the force-elongation measurements done in different biopolymers us- 
ing Atomic Force Microscopes and Optical Tweezers are the best known examples |l4|. This 
new experimental capability has proved that our most advanced models are still incomplete 
and cannot describe modern experimental data adequately. One example of this statement 
is the force-elongation curve of the polyprotein made of the 127 module of human cardiac 

n 

titin 15 . This curve shows discontinuities that occur when one module of the polypro- 
tein unravels. The data is generally fitted to the WCM. But, this model cannot describe 
structures like a-helices or /3-sheets present in polyproteins. 

The solution of the WCM is mathematically involved and the addition of secondary struc- 
ture like a-helices does not simplify the model. On the contrary, the mathematical complex- 
ity of this "extended" WCM is expected to be even more involved. As a consequence of this, 
we have started to explore new mathematical tools that might lead to alternative approaches 
to the models of polymers. One tool that has proved to be potentially useful is the General- 
ized Borel Transform (GBT) which we present in this paper. This technique was borrowed 
from Quantum Mechanics (QM) and Quantum Field Theories (QFT) where it has been 
successfully applied to solve different problems such as heavy quark-antiquark interaction 
in QFT [la] and singular scattering potentials in QM In addition, this mathematical 
tool allows one to obtain analytical solutions of parametric integrals like Laplace/Mellin 
transforms for all range of its parameters. Therefore, it is extremely useful to study non- 
perturbative regimes. In this paper, we show how to use this computational tool in the field 
of statistical mechanics of single polymers by computing exactly the polymer propagator of 
the RFM of linear macromolecules. In the companion paper, we generalize these ideas to 
the statistical mechanics of linear semiflexible polymer chains, as described by the WCM, 
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and present new results for this model. 

The RFM is the simplest attempt to describe the statistical properties of flexible polymer 
chains. Thus, it is the first model where the GBT should be applied. In this model, polymers 
are assumed to be made of n segments of length I (=Kuhn length) which are connected in 
a sequential manner. These segments are free to rotate around any axis as long as the move 
does not break the polymer chain. Figure [T] depicts the RFM. 

This paper is organized as follows. In section II we show how is possible to calculate 
Fourier sin/cos transforms using the Generalized Borel Transform and we present a brief 
summary of the mathematical aspects of this technique. Then in section III we solve exactly 
a particular Fourier sin transform, namely the polymer propagator (normalized distribution 
function) of the RFM. Finally, in section IV we present the conclusions and discuss some 
other models in polymer physics where this new application of the GBT might lead to exact 
solutions. 

II. THE GENERALIZED BOREL TRANSFORM 

In the previous section we have shown that the Fourier Sine transform plays a very impor- 
tant role in the evaluation of the statistical properties of models for single macromolecules. 
Therefore, let us start the description of the GBT by writing the general expression of a 
Fourier Sine transform of a function H (k, a) 

POO 

Q (R, a)= sin (Rk) H (k, a) dk. (4) 
Jo 

Furthermore, consider the Laplace transform of the same function which is given by the 
definition 

POO 

S (g,a) — H (x, a) exp (—gx) dx g > 0. (5) 
Jo 

Then, we observe that we can obtain the Fourier Sine transform, Q (R, a), from the Laplace 
transfom, S (g, a), as the analytic continuation of S (g, a) to the complex plane as follows 

Q (R, a) = Im {S (g = -iR, a)} . (6) 
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Consequently, we will focus on the evaluation of Laplace transforms. For this purpose, we 
will employ the GBT technique described hereafter. 

The main goal of the GBT is to obtain analytical solutions of parametric integrals of 



the Mellin/Laplace typejl^, 12, 3] for all the range of values of the parameters. Therefore, 
this technique is very useful to study non-perturbative regimes. The basic idea of method 
consists of introducing two auxiliary functions, S(g,a,n) and B\ (s, a, n)(the Generalized 
Borel Transform). These functions depend on auxiliary parameters called n and A. These 
parameters have no physical meaning and are introduced for the sole purpose of helping in 
the computation of an explicit mathematical expression for S (g, a) for all values of the true 
parameters g and a. 

Let us start with the mathematical definition of S (g,a,n) which is the following 

/"OO 

S(g,a,n)= / x n H (x, a) exp (— gx) dx. g>0 (7) 
Jo 

We have explicitly extracted a factor x n from the function to be transformed. This integral 
is related to the Laplace transform, Eq. by the following relationship 



S (g,a,n) 



Qn 

dg n 



S (g,a) 



(8) 



which can be inverted to give 



S (g,a) 



n-l 



/ dg--- dg S (g,a,n) + ^2c p (a, n) g p . 

i — , P =o 



(9) 



The finite sum comes from the indefinite integrations. Note that all the coefficients vanish 
whenever the Laplace transform, Eq. fulfills the following asymptotic condition 



lim S (g, a) = 0. 



(10) 



In addition, the expression given by Eq. is valid for any value of the parameter n, 
in particular for n ^> 1 where the GBT provides an approximate analytical expression for 
S (g, a, n) as we describe below. 

Let us define the Generalized Borel Transform of S (g, a, n) as follows 



POD 

B x (s, a,n) = - / exp [s/ r] {g)\ 
Jo 



n -As 



(9) 



+ 1 



S (g,a,n) dg (g) 
[V (9)} 2 dg 



dg, Re(s)<0 (11) 
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where A is any real, positive and non-zero number, and r\ is defined as follows l/rj = 
A(exp(g/A) — 1). Then, B\(s,a,n) is an analytic function for real values of s less than 
zero. Moreover, the analytic continuation to the other half of the complex plane generates 
an analytic function with a cut on the positive real axis. 

In order to invert the transform defined by Eq. ([TT| . we note that the change of variables 
u(g) = l/rj — Aln[l + I/A77] transforms the integral, Eq. (fPT| . into a Laplace transform 

POO 

Bx(s,a,n)= / exp [su] L\ (S, a, n, u) du, Re(s)<0 (12) 
Jo 

where L\(S,a,n,u) depends on S(g,a,n). Consequently, the inverse Laplace transform of 
Eq. (fT2l provides a procedure for the evaluation of S (g,a,n) by integrating B\ (s,a,n) on 
the imaginary axis or over the discontinuity of B\(s,a,n) on the cut. After a change of 
variables we can write S (g, a, n) as follows 

/oo /»oo 
/ exp [G (u>, t, g, A, a, n)] dwdt. (13) 
-00 J — 00 I I 

G (w, t, g, A, a, n) is given by the following expression (for more details see Ref. [16]) 

G (w, t, g, A, a, n) = -s (t) u {g) + 1 - In {r [A (s (t) + x (w))]} 
+ {A [s (t) + x (w)] - 1} In (As (t)) - As (t) + w + In [x (w) n H (x (w))], 
where s (t) — A exp (t) and x (w) = exp (w). 

Note that the expression given by Eq. (JTTfl) is valid for any non-zero, real and positive 
value of the parameter A. However, the resulting expression for S (g, a, n) does not depend 
on A explicitly. Thus, each value of the parameter A defines a particular Borel transform. 
Consequently, we can choose the value of this parameter in such a way that it allows us to 
solve Eq. (13) . 

The dominant contribution to the double integral is obtained using steepest descent! 
in the variables t and w. For this purpose, we first compute the expressions of the saddle 
point t a (g, a, n) and w Q (g, a, n) in the limit A 3> l.The results are the following 



t (g, a, n) = In 



x 2 (g,a } n) 



w (g, a, n) = In [x Q (g, a, n)] , (15) 



/ (x (g,a,n),a,n) 

where x Q (g, a, n) is the real and positive solution of the implicit equation coming from the 
extremes of the function G (w,t, g, A,a, n) in the asymptotic limit in A. Explicitly, the 
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equation is 



x la 2 = f (x , a, n) [f (x , a, n) + 1] , (16) 

where 

/ (x ,o,n) = 1 + ra + x . (17) 

QUiXq 
(the Hessian of G (w, t, g, A, a, n) at the 

saddle point must be positive). Let us call the Hessian D( ). Its mathematical 

expression is 



D (x Q , a, n) = -x ^— ^ — [1/2 + / (x , a, n)} + / (x Q , a, n) [1 + / (x G , a, n)] . (18) 

Observe that in the range of the parameters where / (x Q , a, n) 3> 1, which is fulfilled when 
n > 1, we can keep the second order term in the expansion of G (w, t, g, A, a, n) around the 
saddle point. Then, we can approximate the double integral in Eq. ([T3|) as follows 

n / \ (1 — eX P (~ fi'A)) r^y / > \1 /, ^\ 

^Aprox [9, a, n) = 4tt exp [G (w OJ t OJ g, A, a, ra)J . (19) 

D [x , a, n\ 



In the limit A —>■ oo we obtain the following approximate expression for S (g,a, n) 

Sa P tox (jg, a, n) = V27Te~ 1/2 Q ' = [x ] n+1 H [x a , a] exp [-/ [x a , a, n]] . (20) 

y/D [x , a, n\ 

Note that the expression given by Eq. ([20"|) is valid for functions H (x, a) that fulfill the 
following general conditions. Firstly, the relationship given by Eq. (fT6|) must be biunivocal. 
Secondly, D (x ,a,n) must be positive at x . Thirdly, / (x ,a,n) must be larger than one. 
In particular, this condition is fulfilled when n> 1. These conditions provide the range of 
values of the parameters where the approximate solution, Eq. (j2Hj) , is valid. 

Finally, we replace the expression given by Eq. (ffil into Eq. (jOJ to obtain an approximate 
analytical expression for the Laplace transform S (g,a). In particular, in the limit n — > oo, 
we obtain the following analytical solution for S (g, a) 

lim f f 

S(g,a)= (-) dg - - dg S Aprox (g,a,n) . (21) 

77, — * OO J J 
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One particular case of this result is the one where H (x, a) does not contribute to the 
saddle point. This is the case when / (x , a, n) can be approximated by 1 + n (the derivative 
of H (x Q , a) is negligible). Then, the saddle point solution is x Q (g, a, n) ~ (n + 3/2) fg and 
the expression of SUprox (flS a i n ) is 

r (n + 1) 

SAprox (g,a,n) ~ — — H [x (g, a, n) , a] n > 1. (22) 

Another important property of the expression given by Eq. (|2T| is that, in the limit 
n — >• oo, the approximate solution, Eq. ([20]). becomes an exact solution for Eq. 0. Thus, 
as long as the n indefinite integrals are calculated without approximations, then Eq. (|2*T|) is 
an exact solution for Eq. JSJ). 

In summary, the procedure to use the GBT to compute Fourier Sine transforms is the 
following. Firstly, one has to solve the implicit equation, Eq. (fT6|) . for n > 1 to obtain 
the mathematical expression of x (g,a,n). Replacing this expression into Eq. (|2T)|) and 
doing the n indefinite integrals in Eq. (|2T| . we get the exact/approximate expression for 
S(g,a). Finally, one has to compute the analytic continuation of S(g,a), Eq. JHJ), to get 
the exact/approximate solution of the Fourier Sine transform, Eq. (J3J. 

In the next section we apply this technique to solve exactly the statistical mechanics of 
flexible macromolecules. 



III. APPLICATION 

Let us start by analyzing the polymer propagator predicted by the Random-Flight Model 
which is given by 

/n / n \ 

d {K k } J] r (Rj) 8 IJ2 R i - R ' ( 23 ) 
i=i \i=i / 

where Rj is the bond vector between the (j — l)-th and j-th beads, n is the total number of 
segments per polymer chain, R is the end-to-end vector and r (Rj) is given by the formula 

T{Rj) - 47TZ* ■ (24) 
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Using the Fourier representation of the delta function, we obtain 



P (R, n) 



d 3 kexp (—iH ■ k) 
(2tt) 3 (AnP) n 
d 3 k 



/n / n \ 

rf{R fc }n 5 (l R J |- / )exp^^R J -k 
3=1 \ j=l J 



(25) 



(27T)'- 



exp (— zR • k) .fT (k, n, /) 



where the characteristic function, K (k,n,l), is 



K(k,n,l) 



sin(|k|/) 
Ikl/ 



(26) 



The evaluation of the angular integrals in Eq. |2*K|) is straightforward. After rescaling R 
and k with the Kuhn length, /, we obtain the final expression for the polymer propagator 



P{R,n) 



dk 



sin(^)(^l k 



(27) 



where k = |k| and R = |R|. 

This integral representation of the polymer propagator is a Fourier Sine transform and 
can be solve exactly using GBT. Then, our first step consists of expressing the polymer 
propagator, Eq. (|27|) . in terms of a Laplace Transform. For this purpose we define the 
function 



S (b, n) = 



dw 



exp (—wb) I w 

\ w 



(28) 



from where we recover the expression of the polymer propagator, Eq. (|27j) . as the analytic 
continuation of the function S (6, n) to the complex plane 



P (R, n) 



o — I m \S (b = —iR, n)} 

(2n) 2 R 1 V n 



(29) 



Let us now rewrite Eq. ()28|) as follows 



Qn 

S(b,n) = — 

V > J Q c n 



dw 



wexp (—wb) exp I c 



sin l w 



w 



d" 



c=0 



c=0 ' 



(30) 



where 



GA (b, c) 



dww exp (—wb) H (w, c) 



(31) 
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and 



sin (w) 



w 



(32) 



H (w, c) = exp 

Then, the integral expressed by Eq. (pHT) satisfies all the requirements of the GBT technique 
Consequently, we evaluate this integral in the following way 



GA (b, c) 



lim 

iV -> oo 



db / db GA 



N 



where we have defined 



(33) 



GA N {b, c)= dw [w N+1 exp {-wb) H (w, c)] . (34) 
Jo 

In the asymptotic limit of N —>■ oo, the GBT provides an analytical solution for Eq. (|3~4T) . 
Following the technique, we first solve the following implicit equation for w , Eq. (|16|) . 



N + 1 + Wo ^- In [H (w , c)] | S.N + 2 + w ^- In [H (w , c)} \ = w 2 a b 2 , (35) 



whose asymptotic solution is 



N + 5/2 Ar 1 
w ~ — iV > 1. 



(36) 



Replacing this expression for w Q in the expression provided by the GBT, Eq. (J2*0*J) . we obtain 



^ft^IM.p.) iV » 1 . ,3T) 

Furthermore, we replace Eq. (I37}l into Eq. (j3~31) and the resulting expression into Eq. ([371 . 
then we obtain 



lim * lim i-f db U^H^c) 

S(b,n)= —{ I . }. (38) 

We now proceed to exchange the order of the operators, first we evaluate the n-th deriva- 
tive of the function H with respect to c and, afterward, we take the limit of c — > . As a 
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result, we obtain 



Next, we solve the N integrations. Using properties of the function sin (x) we can write 
S (b, n) for any odd number of segments as follows 

S (b, n) = -L J2 (") V+fc \]M{N,n, k, b) , (40) 



2 

k=0 



k 



where 



M(N,n,k,b) = ^ H^E <i«^ / * . . ./ (41) 



We note that the only powers of b in Eq. (jUj) that fulfill the asymptotic behavior of the 
function S(b,n), Eq. (fTOjl. are those that satisfy the condition r > (n — 1). Consequently, 
the N indefinite integrations are exactly doable, the result is 

f 1 T(2 + r-n) (-) N 

db " " " J db b N+2-n+r = T ( N + 2 - n + r) b 2+r ~ n ' ^ 



Replacing Eq. (|42j) into Eq. (|4Tj) and introducing the dummy variable r = x + n— lwe can 
write 



M ( w^(^- 2 *)r'x;(i^)^ 



lim N^T (N + 2) 

X AT^oo T(N + x + l) ' 

which, after using the asymptotic properties of the Gamma function [3] 

lim N^T (N + 2) 

iv^oo r(iv + x + i) " ' 

becomes 



(43) 



(44) 



M („, ifc, 6) = \lm £ (i (n - 2k))^ ' £|±l) . (45) 



11 



The sum over x is doable, the result gives the following expression for M (n, k, b) 



M (n, k, b) = -Im [(i (n - 2k)) n ~ l FD (n, k, b)] , (46) 



where we have defined 



FD (n, k, b) = ^ 



(n-2k) 



v^r (n) / 47 ^ 

-2^)3^2 ([1,1, |] ,[^,^] 



i[n — Ik) a 2 ' ' 2J ' L 2 ' 2 J ' &2 
+ 6 r(n + l) 

3 F 2 ([,,],[,] ,x) is the Generalized Hypergeometric function [20]. From Eq. ijlfij) we can see 
that the imaginary part affects only the function FD (n,k,b). Thus, we obtain the final 
expression for S (b, n) 



" ( "' ^ = £ H [k J + 1) ' (48) 

The last step to obtain the analytical expression of the polymer propagator consists of 
inserting Eq. P%Jl into Eq. (j2*9*j) and computing the analytic continuation of the resulting 
expression to the complex plane through the substitution b = —iR. After doing these 
computations, we arrived at the following expression for the polymer propagator 



mn) = — ^>>r +J I . \{n-2kf 1 



xlm 




r (n + 1) 

(49) 

n-2ky 



R 2 



2 ' 2 

This expression can be simplified even further if we use the well known analytical prop- 
erties of the Hypergeometric function 3 F 2 (z) |2^| which is an analytic function for values of 
\z\ < 1 and its continuation to the rest of complex plane generates one cut on the positive 
real axis starting at Re (z) = 1. This implies that only values of ^ n ~^ > 1 will contribute 
to the imaginary part of 3F2 (z). Consequently, this condition reduces the number of terms 
in the fc-sum such that the last term of Eq. fi9j) is k = ■ 

The explicit evaluation of Im | 3 F 2 (Jl, 1, §] , p±i, ^f] , ^=^) } can be found in Ap- 
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pendix A. The final expression is 



Im 



3-^2 



1,1,?' 



R 2 ir 



n+1 n + 2 



2 ' 
T(n + 1) 



[n-2k- R] 




n>2 



(50) 



2{n-2k) n Y{n-l) 

Finally, we replace Eq. (f50|) into Eq. fi9|) to obtain the exact expression for the polymer 
propagator 



P (R, n) 



1 



2 n+1 nR 



£ 

fc=0 



n 



k 



[n-2k- R] n ~ 2 
T(n-l) 



(51) 



Equation (|5Tj) is valid for odd number of segments but, it is extended to polymers with 
any number of segments larger than two via ana vtic continuation. 

Therefore, we have obtained the well-known 2l|, |22| exact analytical expression for the 
polymer propagator of flexible chains, Eq. lj27j) . with any number of segments, n, and any 
end-to-end distance, R. 

Observe that the result given by Eq. (JBTj) can be used to describe the statistical properties 
of polymers with other topologies. For example, consider the case of a flexible m-arm 
star polymer as shown in Fig. . Since the polymer is flexible, then each arm behaves 
independently from the other ones except for the fact that all of them start at the origin. 
Thus, the probability of finding the end of the j-th arm in the shell of radius Rj with 
thickness dRj centered at the origin is 



AnR]P (Rj.ji^dRj, 



(52) 



where rij is the number of segments in the j-th arm. If we consider all the arms, then the 
probability of finding the end of the first, second, etc. arms in the shells of radii Ri,R2, etc. 
with thicknesses dR\, gLR 2 , etc. centered at the origin is 



(47r) m JJ J RjP( J R i ,n i )d J R i . (53) 

Other probability distributions for star polymers can also be computed easily. 

Another example is the case of ring (cyclic) polymers. Figure shows this topology. From 
this figure and following the steps presented in this paper for linear polymers it can be 
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proved that the probability of finding any pair of segments separated by a distance R should 
be proportional to the product of two propagators of the form given by Eq. I)51jl 

P Rmg (R, s,L-s)^P {R, s) P{R,n-s), (54) 

where n is the total number of segments in the ring and s is the number of segments (along 
the contour of the polymer chain) between the two chosen segments. 

The aformentioned two examples clearly show that the results obtained for linear poly- 
mers using the GBT can be used for polymers with other topologies thus, increasing the 
range of models that are mathematically tractable with the GBT. 

IV. CONCLUSIONS 

In this paper we have described a new mathematical method called the Generalized 
Borel Transform and applied it to compute some statistical properties (polymer propagator) 
of models of flexible polymers. Specifically, we showed how the GBT is constructed and how 
to use it to compute Mellin/Laplace transforms. Moreover, some mathematical properties 
were presented. The application of this technique to the statistical mechanics of single 
flexible polymers led to the exact solution for the polymer propagator of linear polymers. 
The propagator obtained turned out to be a finite sum of polynomials valid for any end-to- 
end distance, R, and number of segments, n. Furthermore, this result was used to construct 
distribution functions for two other topologies, rings and stars. 

The exact computation of the polymer propagator of the RFM is a straightforward cal- 
culation that requires simple mathematics when the GBT is used. This mathematical sim- 
plicity of the GBT makes it a potentially very useful computational tool for more complex 
models of single polymer chains because it does not add any complexity to the physics of 
the starting model. 

Equation (|5T| together with its extensions to stars and rings, Eqs. (|53|) and (f5"i|) . and the 
discussion presented in the introduction show that the GBT can solve exactly a wide range 
of models for polymers. However, more advanced models of single polymer chains like the 
Wormlike Chain Model or helical polymers where the bond vectors are correlated with each 
other through potential interactions are not exactly soluble with the GBT at present. A 
generalization of the GBT to multidimensional integrals is required to address these models. 
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Appendix A: EVALUATION OF THE IMAGINARY PART OF THE HYPERGE- 
OMETRIC FUNCTION 



In this appendix, we calculate the expression Jm{ 3 F 2 (z)}. For this purpose, we use the 
following integral representation of the Hypergeometric function 



3-^2 



— V 



A A+ 1 



A + /i A + /i + l 



-1 



(Al) 



x f Q x^ 1 [q - xf' 1 \x 2 + w 2 ] v dx, A, a > 0, Re (-) > 
Jo |— | W 

where B (A,//) is the Beta function |2f|. 

We now assign the values v = — 1, A = 2, /i = n — 1, q = n — 2k, w = b to the parameters 
in Eq. (|A1|) to obtain 



3-^2 



1,1,- 

' '2 



n + 1 n + 2 



2 ' 2 

n-2fc 



(n - 2fc) 2 



x 



(A2) 



o 



b 2 j b- 2 (n-2k) n B(2,n- 1) 
x [n — 2k — x] n ~ 2 [x 2 + 6 2 ] 1 cfx. 

This integral representation is valid only for n > 2. Therefore, when we take the analytic 
continuation to the complex plane as before {b = —iR), we can express the imaginary part 
of the Hypergeometric function as follows 



Im < 3 F 2 



1,1, . 



n+1 n+2 
2 ' 2 

n-2k 



(n - 2fc) 2 



i? 2 



i? 2 



[n-2k) n B (2,n- 1) 



xlm 



x In — 2k — x 



n-2 r 2 







[x 2 -/? 2 ]^. 



Thus, we need to evaluate 



rn— 2k 

L = Im / x [n — 2/c — x] n ~ 2 [x — R]~ [x + R]~ x dx. 
Jo 



(A3) 



(A4) 
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After analyzing the analytical behavior of the integrand, we concluded that we can exchange 
the operations of integration and imaginary part to obtain 



n-2k 



x [n — 2k — x} n 2 [x + R] 1 Im { [x — R] 1 } dx. 



(A5) 



Thus, we have to compute 



LS = Im 



1 



(A6) 



{x-R) 

first and, afterward, we have to solve the integral given by Eq. (IA5J) . 

The analytical behavior of the function (x — R)~ is well known. It is an analytic function 
for \x\ > R but, its analytic continuation to the complex plane generates a cut on the real 
axis in the range —R < Re (x) < R. This cut generates its imaginary part which is 



Im 



1 



ttS (x — R) . 



(A7) 



(x-R) 

Thus, replacing the expression given by Eq. (|A7|) into Eq. (|A5|) and making the change of 
variables y = x — R, we obtain 



where 



rn-2k-R 

L = ~ / F k (y, n, R) 5 (y) dy, 

-R 



F k (y, n,R) = {y + R)[n-2k-R- y] n ~ 2 [y + 2R] 1 



The result of the integration gives 



(A8) 



(A9) 



n-2 



L = - [n - 2k - R] 
Finally, we replace Eq. (|A10J) into Eq. (|A3|) to obtain the final expression 



(A10) 



Im < 3 F 2 



1,1, 



n + 1 n + 2 



(n - 2kf 
W 2 



TT 



R? 



2 (n-2k) n B(2,n-l] 



[n-2k- R] 



n-2 



(All) 
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List of Figures 

FIG. EJ Linear, ring and 4-arm star topologies for flexible polymers. R indicates the end-to- 
end vector for the linear topology and the relative position of two segments in the case 
of a ring polymer. Rj indicates the position of the end of the j-th arm in the star 
topology. 

FIG. [2 The Random Flight Model of polymer chains. / is the Kuhn length and 9 is the bond 
angle. 
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Linear Ring 4-arm Star 

FIG. 1: Linear, ring and 4-arm star topologies for flexible polymers. R indicates the end-to-end 
vector for the linear topology and the relative position of two segments in the case of a ring polymer. 
Hj indicates the position of the end of the j-th arm in the star topology 
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FIG. 2: The Random Flight Model of polymer chains. I is the Kuhn length and 9 is the bond angle. 
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